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ABSTRACT: This paper proposes a framework for the analysis and design optimization of models subject to 
parametric uncertainty where design requirements in the form of inequality constraints are present. Emphasis is 
given to uncertainty models prescribed by norm bounded perturbations from a nominal parameter value and by 
sets of componentwise bounded uncertain variables. These models, which often arise in engineering problems, 
allow for a sharp mathematical manipulation. Constraints can be implemented in the hard sense, i.e., constraints 
must be satisfied for all parameter realizations in the uncertainty model, and in the soft sense, i.e., constraints 
can be violated by some realizations of the uncertain parameter. In regard to hard constraints, this methodology 
allows (i) to determine if a hard constraint can be satisfied for a given uncertainty model and constraint structure, 
(ii) to generate conclusive, formally verifiable reliability assessments that allow for unprejudiced comparisons 
of competing design alternatives and (iii) to identify the critical combination of uncertain parameters leading to 
constraint violations. In regard to soft constraints, the methodology allows the designer (i) to use probabilistic 
uncertainty models, (ii) to calculate upper bounds to the probability of constraint violation, and (iii) to efficiently 
estimate failure probabilities via a hybrid method. This method integrates the upper bounds, for which closed 
form expressions are derived, along with conditional sampling. In addition, an formulation for the efficient 
manipulation of hyper-rectangular sets is also proposed. 


1 INTRODUCTION 

Design under uncertainty arises in numerous disci- 
plines including engineering, economics, finance and 
management. Achieving balance between robustness 
and performance is one of the fundamental challenges 
faced by scientists and engineers. 

Literature in probabilistic control (Crespo and 
Kenny 2005), stochastic programming (Kail and 
Wallace 1994) and stochastic approximations (Er- 
moliev 1983) provide several mathematical tools for 
Reliability-Based Design Optimization (RBDO). The 
algorithms at disposal can be classified according to 
the way they integrate inequality constraints. Hard 
constraints (Darlington et al. 1999) are those that 
must be satisfied for all possible realizations of the 
uncertain parameter. Strategies to solve the resulting 
semi-infinite optimization problem usually require 
nested searches for the identification of the worst- 
parameter-realization (Darlington et al. 1999). On the 


other hand, Soft constraints (Acevedo and Pistikopou- 
los 1998; Samsatli et al. 1998; Rustem and Nguyen 
1998; Royset et al. 2001 ; Crespo and Kenny 2005) are 
those that are not required to be satisfied for all pos- 
sible realizations of the uncertain parameter. They re- 
sult from slack design requirements, overly large un- 
certainty sets or limiting design structures. Chance- 
constrained programming (Kail and Wallace 1994), 
sampling-based techniques, asymptotic approxima- 
tions (Rackwitz 2001) and penalty-based optimiza- 
tion (Rustem and Nguyen 1998; Kail and Wallace 
1994) are some of the strategies most commonly used 
to tackle this problem. 


This paper proposes strategies for robustness anal- 
ysis and RBDO of problems subject to hard and/or 
soft constraints. A unifying mathematical framework 
is first introduced. 
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2 FRAMEWORK 

Definition 1 (Uncertainty Model). If p G M dim ( p ) is 

an uncertain parameter, its uncertainty model is given 
by the support A p C M' ]mi!p) , and a designated point 
P e A p . 

The uncertainty model is specified by the ana- 
lyst/designer. The intent is that the support set A p be 
chosen so that the actual value of the parameter p lies 
somewhere within it. A realization of the uncertain 
parameter is a value of the parameter selected from 
the support set. The designated point, whose selec- 
tion is subjective, can be interpreted as the realization 
that best represents p. 

Consider now the situation that a system depends 
on the uncertain parameter p G and the de- 

sign variable d G M dim(d) . Suppose that g : M < ' lllTl(p - 1 x 
]g>dim(d) ]g>dim(g) j s a set 0 f cons t ra int functions on 
the system (with the convention that positive values 
represent constraint violations). If these are consid- 
ered hard constraints, the system corresponding to 
given values of d and A p will be judged acceptable 
if g(p, d) < 0 \/p G Ap. 

For a given vector g of constraints and a given 
design d, denote the failure (or constraint violation) 
set by F(d,g) = (j d ™ (9) Ffid,g) where Tfid.g) = 
{p '■ 9 i{P>d) > 0} is the set of p- values where con- 
straint number i is violated. In a context where g and 
d are understood, lF(d,g) and SFfid, g) will be rep- 
resented simply as T and T t . Using the symbol d to 
represent the topological boundary of a set, if the con- 
straint functions are continuous, each p G DT, satis- 
fies gfip. d) = 0 and eachp G dT satisfies g j(p, d) = 
0 for some j 1 . 

The manifold dT commonly referred as the limit 
state surface, separates the parameter space into two 
regions: one where all the constraints are satisfied and 
one where they are not. For a fixed design d, the fail- 
ure region T either overlaps the support set A p , case 
in which the design d is called non-robusf, or T and 
Ap are disjoint, case in which the design d is called 
robust. In the latter case, the degree of robustness can 
be quantified by measuring the separation between the 
two sets. 

Selecting the support set usually involves some en- 
gineering judgment. One reasonable choice might be 
to confine each component of p to a bounded interval. 
This leads to the choice of A p as a hyper-rectangle. A 
natural choice for the designated point p is the geo- 
metric center of A p . If m is the vector of half-lengths 
of the sides of the hyper-rectangle, the rectangle is 

'The converse is not universally true. We will assume in this 
paper that the boundary of T t is exactly described by the points 
where that constraint function assumes the value of zero; i.e., 
dTj — {p : g t {p, cl) = 0}. We believe that in a realistic engi- 
neering setting, the loss of generality in making this assumption 
is unimportant. 


represented by the notation 7 Z p (p, m ), and defined by 
n p (p,m) = {p : p i G [pi - m i ,p i + mj, 1 <i< 
diiii(p ) }. If vector inequalities (a < b) and vector ab- 
solute value (|a|) are understood to hold component- 
wise, this hyper-rectangle is given by {p : \p — p\ < 
m} where m > 0. 

Another reasonable choice for the support set is a 
hyper-sphere. The hyper-sphere of radius R centered 
at p will be denoted as S p (p, R ). 

Definition 2 (Homothetic Sets). Two sets A and B 
are homothetic with respect to the homothetic center 
x at a similitude ratio of a if B = {b : b = a (a — x) + 
x, a G A}. Such a B is completely determined by A, x, 
and a. We use the notation B = H(A, x, a) to express 
this relationship. 

B = Ti(A, x, a) means that B can be created from 
A by forming every vector from the fixed point x to 
each point of A; stretching or shrinking it by a factor 
of a; and, with the roots of the stretched or shrunk 
vectors all fixed at x, collecting all the new endpoints 
to form the set B. 

For purposes of this paper, two uncertainty mod- 
els will be called proportional if they have the same 
designated point and they are homothetic with re- 
spect to that designated point as homothetic center. 
This means that one of the two support sets can be 
formed from the other by expansion or contraction 
by some positive factor, which is the similitude ratio, 
about the common designated point. For instance, the 
hyper-rectangles 7 Z p (p, m ) and 7Z p (p. am) are pro- 
portional sets for a > 0. 

The notions of Critical Parameter Value (CPV) 
and Parametric Safety Margin (PSM) are now intro- 
duced. For simplicity sake, the presentation of the ma- 
terial will concentrate on the case where the desig- 
nated point is in the non-failure region. Intuitively, 
one imagines that a set proportional to the support 
set of the uncertainty model is being expanded ho- 
mothetically with respect to its designated point until 
its boundary just touches the failure region. The point 
where the expanding set touches the infeasible region 
is the CPV. The PSM is a metric that quantifies the 
size of the set proportional to the support set that has 
the CPV on its surface. The mathematical background 
for these notions is presented next. 

Definition 3 (Critical Similitude Ratio). Let d be a 

fixed design and let A p and p be the support set and 
designated point of the uncertainty model of p. Sup- 
pose that g : M dim ( p ) x M dim ( d ) -> R d Ma) i s a set of 
constraint functions. The set of expansion/contraction 
factors which produce homothets of A p with respect 
to p which contain constraint violation points is de- 
noted 

A(p,A p ,d,g) = ja > 0 : Tt p (A p ,p,a) f 0j . 


2 



The greatest lower bound ofA(p, A p , d,g), given by 
a(p,A p ,d,g) = inf ^A(p,A p ,d,g)^j (1) 

will be called the Critical Similitude Ratio. 

Arguments that can be understood from context 
will be omitted. 

Definition 4 (Critical Parameter Value). Let d be 

a given design and let A p , p, and g define the un- 
certainty model and the set of inequality constraints. 
Any p G 'H P (A P . p. a) on a limit state function, i.e. 
gfp, d) = 0 for some i, will be called the Critical Pa- 
rameter Value. 

Some properties of the CPV are provided in ref- 
erence (Crespo, Giesy, and Kenny 2006). Whenever 
A p is compact and the g are continuous functions, 
the existence of a CPV is guaranteed. Note also that 
the CPV might not be a realization of the uncertain 
parameter, i.e., p might not belong to A p . 

Formal definitions of the PSM for hyper-spherical 
and hyper-rectangular supports are provided in Sec- 
tion 3, as are expressions for the calculation of these 
PSMs. In general, the PSM is a measure of the ro- 
bustness of the design d to uncertainty in p that will 
assume non-negative values in the Feasible Design 
Space (FDS). The FDS is the set of designs satisfying 
all the constraints at the designated point. The PSM 
is proportional to the degree of resilience of d to un- 
certainty in p. If the PSM assumes the value of zero, 
there is no resilience since at least one of the con- 
straints is active for p , i.e. there exists an infinites- 
imally small perturbation from p leading to a con- 
straint violation. Negative PSM values are attained by 
designs outside the FDS. On the other hand, the PSM 
is infinite when g is insensitive to p. In general, the 
PSM is uniquely specified by the design point d, the 
support set A p , the designated point p. the constraint 
functions g, and the critical similitude ratio a. 

Definition 5 (Robust Design Space). The Robust 

Design Space (RDS) is the set of designs dfor which 
g(p, d) < 0 for all realizations of p. Each member of 
this set is called a robust design. 

In this paper, robustness refers to the ability of a 
given design to satisfy the design requirements pre- 
scribed in g < 0 given an uncertain model for p. In 
principle, such requirements can be enforced by ei- 
ther hard and soft constraints. Developments for the 
former type are presented next. 

3 HARD CONSTRAINTS 

Problem Statement: evaluate whether the design d is 
robust for the support set A p and the hard constraint 

g(p d) < o. 


In general, the support set can have any shape 
and size. We assume hereafter that A p is a compact 
set with a designated point p. A general robustness 
test is available in (Crespo et al. 2006). Support sets 
with hyper-spherical and hyper-rectangular support 
geometries are studied next. 


3.1 Hyper-Spheres in p- space 
Hyper-spherical supports result from uncertainty 
models where norm bounded perturbations from the 
nominal parameter value p are allowed. This implies 
that p is the geometric center of the hyper- sphere. 
Problems with this class of support sets are the sim- 
plest since the CPV is calculated by solving a mini- 
mum norm problem in p-space. The CPV for A p = 
S p (p , R) is given by 


p = argmin 

p 


dim(g) 

P-p || : n 9k{p,d) =0 

k= 1 


( 2 ) 


Definition 6 (Spherical PSM). The Spherical Para- 
metric Safety Margin corresponding to the design d 
for A p = Spip. R ) is defined as 

Ps(p,p,d) = ||p-p||, (3) 

where p is the corresponding CPV and p, the des- 
ignated point of the uncertainty domain, satisfies 
g(p, d) < 0. Note that ps = ocR. 

Lemma 1 (P-Test). Let g < 0 describe a set of in- 
equality constraints and A p = S p (p, R) be the sup- 
port of the uncertain parameter p. The design d is 
robust if and only if ps > R. 

The exact boundary of the RDS corresponding to 
a hyper-spherical support set is confined by an iso- 
spherical-PSM manifold, i.e. surface in the design 
space prescribed by p s = c where c is a constant. 


3.2 Hyper-Rectangles in p- space 
When the components of the uncertain parameter are 
independently bounded the support set is a hyper- 
rectangle. Note that this geometry allows for the ma- 
nipulation of uncertain parameters having different fi- 
delities. If p is the geometric center of the set, such 
a hyper-rectangle is 1Z p (p,m). A robustness test for 
this geometry is introduced subsequently. The mathe- 
matical background for this is presented next. 

Definition 7 (Q-Transformation). Let Ap = 

lZ p (p,m) be the support of the uncertain pa- 
rameter p. The Q-Transformation, denoted as 
q = Q(p,p, m), and given by 


Q(P,P,m ) 


A max{|fc|}fc 


(4) 
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where 

k = diagjm}^ 1 (p — p) , 

transforms A p into a unit hyper-sphere in q-space. 
The inverse transformation, p = 0 '(q.p.rn), is 
given by 

Q~\q,P,m) = p + diag{m} . (5) 

max{ \q\\ 

The Q-Transformation maps hyper-rectangles pro- 
portional to 7 Z p (p,m) into hyper-spheres propor- 
tional to <S q (0, 1). The connectivity of the uncer- 
tainty set is preserved by the transformation. Notice 
however, that differentiable functions on p- space are 
transformed into functions on q-space which can have 
derivative discontinuities at points corresponding to 
p- space points where the faces of the homothets meet. 
The notation Q(p) will be used when the arguments 
p and m are understood by context. 

The Q-Transformation allows identification of the 
corresponding CPV by solving the following mini- 
mum norm problem in q-space 


p = argmin 

p 


dim(g) 

Qip) II : n 9k(p,d) 

k= 1 



Definition 8 (Rectangular PSM). The Rectangular 
Parametric Safety Margin corresponding to the de- 
sign dfor A p = 7 Z p (p, m), is defined as 

pn{p,P,m,d)=\\Q{p)\\\\m\\, (7) 

where g(p,d) < 0. 


by discontinuities in the gradient of the Q transforma- 
tion. 

Recall that the infinity norm in a finite dimen- 
sional space is defined as ||£c||°° = sup, : {|x,;|}. Let 
us define the m - scaled infinity norm as ||a?||“ = 
supj{|ccj|/mj}. A distance between the vectors x and 
y can be defined as || cc — y||“. Using this distance, 
the unit ball centered at p is just lZ{p, m). 

The problem of finding, for a fixed design d, the 
CPV for the vector g(p,d ) of constraint functions 
and the uncertainty model with designated point p be- 
comes the problem of finding a vector p of minimal 
distance in this m-scaled infinity norm from p such 
that p touches the failure set. This can be expressed in 
the form of a constrained optimization problem 

p = argmin (||p — p||“ : p E dJ 7 } 

p 

This optimization problem is restated as p = p, where 
i = argmin {HPj.-pII*}, (8) 

3 

P 3 = argmin { ||p - p||“ : P E dJ= 3 ) . (9) 

p 

That is, the CPV problem is solved for each individual 
constraint function, and the answer is selected which 
is closest to the designated point in the m-scaled in- 
finity norm. It should be pointed out that, in the pre- 
ceding discussion, if the m-scaled infinity norm is re- 
placed by the standard Euclidean norm, what results 
is a formulation of the CPV calculation for spherical 
uncertainty domains. 

Assuming that g{p,d ) < 0 and dT { = 
{p '■ 9i(Pi d) = 0}, the problem in Equation (9) 
can be formulated as 


H p (1Z,p , a) is tightly bounded by an hyper-sphere 
of radius equal to p- R centered at p. 

Lemma 2 (Q-Test). Let g < 0 describe a set of in- 
equality constraints, and Ap = 1Z p {p, m) be the sup- 
port of the uncertain parameter. The design d is ro- 
bust if and only if pn > ||m||. 


Pj = argmin {| |p - p||“ : g 3 {p,d) = 0}. 

p 

Rewriting this using the definition of the m-scaled 
infinity norm gives 

• / \Pi~Pi\ f i\ n \ 

p : = argmin < max : gAp.d) =0 > . 

1 p ( l<*<dim(p) 77lj J 


The CPVs in Equations 2 and 6 are p values, not 
necessarily realizations in A p , on the intersection be- 
tween the surface of the homothetic set and the failure 
region. As before, the exact boundary of the RDS for 
hyper-rectangular support sets is prescribed by an iso- 
rectangular-PSM manifold. Relevant aspects related 
to the existence of the RDS are considered in (Cre- 
spo, Giesy, and Kenny 2006). 

3.2.1 Infinity Norm Formulation 
An alternate way to search for the CPV for hyper- 
rectangular supports is presented here. This formu- 
lation allows us to circumvent the problems caused 


The “max” can be eliminated and the objective func- 
tion made differentiable by introducing the similitude 
ratio a defined before 

( Pji&j ) = argmin{a; : q^p, d) = 0, |p* -p*| < am*}, 

p,a 

where 1 < i < dim(p). The non-differentiabilities 
in the objective functions can be eliminated and 
the problem turned into an “inequality constraint 
only” optimization problem which is more “optimizer 
friendly” by changing the constraint on g • from = to 
> (assuming that the constraint functions g are con- 
tinuous) since the optimum must occur on dJ 7 . 
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4 SOFT CONSTRAINTS 

Up to this point, only hard constraints have been con- 
sidered. In practice, there might be cases in which 
the design architecture, i.e., the form in which g de- 
pends on d, or the size of the support set, lead to 
an empty RDS. In such cases, soft constraints allow 
for the search of designs that minimize the severity 
of the constraint violation. Deterministic (Rustem and 
Nguyen 1998) and probabilistic (Royset, Kiureghian, 
and Polak 2001; Crespo and Kenny 2005) approaches 
can be used to tackle problems where partial fea- 
sibility is satisfactory. Non-probabilistic uncertainty 
models have been used thus far. This has been the 
case since the enforcement of hard constraint makes 
the probabilistic information inconsequential. In what 
follows, probabilistic uncertainty models are used. A 
probabilistic uncertainty model is fully prescribed by 
a joint probability density function of p , denoted by 
f p (p), supported in A p . 

Strategies for the efficient solution of RBDO prob- 
lems that involve the estimation of the failure prob- 
ability P[J 7 (d,g)] are studied in this section. Notice 
that if the probability of constraint violation can be 
eliminated, the tools for hard constraints must be used 
instead. The implementation of soft constraints must 
be preceded by determining that the RDS is an empty 
set. Failure to recognize this will lead to ineffective 
numerical implementations. 

By definition, the homothetic sets corresponding 
to the critical similitude ratio are fully contained in 
the non-failure region. If a probabilistic uncertainty 
model for p is available, the probability of being 
within such a set can be quantified. Obviously, this 
probability depends, among several other factors, on 
the geometry of the support set. In order to maximize 
the probability of being within such a homothet, ex- 
pansions/contractions in the standard normal space u 
are desirable. The reader can refer to (Rackwitz 2001) 
for a review on the transformation from p-spacc to u- 
space. We will use u = U (p) to refer to the probabil- 
ity preserving transformation. Most of the concepts 
introduced earlier can be naturally extended to this 
setting with one notable exception. Since the trans- 
formation of the support set A p covers the entire u- 
space, the concept of proportionality between the ho- 
mothets and the support set is no longer attainable. 
This gives us the freedom to select dilating sets of ar- 
bitrary geometry in rr-space. As before, we will con- 
centrate on hyper-spherical and hyper-rectangular ge- 
ometries. We will denote with u the designated point 
in standard normal space. Analogous to the spheri- 
cal and rectangular PSMs concepts used for hard con- 
straints, are the spherical- and rectangular-Reliability 
Indices (RI) for soft constraints. 


4.1 Hyper-Spheres in u-space 
The CPV for this type of dilating sets is calculated 
by solving a minimum norm problem in u- space. In 
particular, the CPV is given by 


{ dim(g) 

||w — w|| : JJ g fU^ 1 (u) , d) = 0 
*=1 

Definition 9 (Spherical RI). The reliability index 
corresponding to the design d for dilating spheres 
centered at u, called the Spherical Reliability Index, 
is defined as 



Ps(u,u,d) = \\ii-u\\, (10) 

where g(u , d) < 0. 

If u = 0, Ps is the conventional reliability index 
and the CPV is the Most Probable Point. With this 
information at hand, a bound to the failure probability 
can be derived. 

Theorem 1 (Bound to P[P\: Hyper-Spheres). The 

number v.>s, given by 


'fis = 1 - erf(y) + \j -e 7 /j (dim (u),p s ), (11) 


' n - 3 
2 


1 n — 1 
2 


fi(n,r) = ( Y[h(j) ] 
J =1 


k= 1 


l -~k 


„2k— 


1 ]f h{i) 


i = 1 


when dim ( u) is an odd number, and by 

'fis = f 2 (dim (u),Ps) , (12) 


/f-1 > 

, 1 n/2 

f- fc 

n h u) 

V =1 ) 

E 

' k= 1 

r 2 M H h(i) 

i= 1 


when dirn(it) is an even number, is an upper bound 
to the failure probability P[T ] corresponding to a 
hyper- spherical set with designated point u = 0 and 
spherical reliability index Ps- In these expressions, 
h(i ) = n — 2i and y = Ps/ v/2- 

4.2 Hyper-Rectangles in u-space 
The developments that follow are an extension of the 
ideas presented in Section 3.2. The CPV resulting 
from dilating a hyper-rectangle centered at u having 
the half-length vector m is given by 


u = argmm 

U 


dim (g) 

\u-u\\Z- n 9i(U~ 1 (u),d)=0 


5 



An equivalent formulation based on the Q- 
Transformation requires Q(u,u,m ) as the 

cost function. In this context, q = Q(u,u,m) 
maps hyper-rectangles in r/, -space proportional to 
7 Z u (u, m ) into hyper-spheres in q-spacc proportional 
to <S g (0,l). For the reminder of this section, Q(u ) 
refers to Q(u, u, m). 

Definition 10 (Rectangular RI). The reliability in- 
dex corresponding to the design d for a dilating 
hyper-rectangle with designated point u and the half- 
lengths vector m, called the Rectangular Reliability 
Index, is defined as 


Pn(u,u,d) = 


Q 


-i 


m i 


\Ali 


— u 


im u 


(13) 


where g(u, d) < 0 and q = Q(u). 

Theorem 2 (Bound to P[JF]: Hyper- Rectangles). 

The number V'-r,, given by 


dim (u) 

f’n — [1 + -’ (Ji) - $(ttj + af)] (14) 

2=1 

is an upper bound to the failure probability P[P\ for 
a hyper-rectangular set with designated point u and a 
rectangular reliability index (3 k- In this expression, $ 
is the cumulative distribution function of a standard 
normal variable and cr, = 6-f’rn l /\\rn\\. 

Note that as long as the search for the CPV con- 
verges to the global minimum, the upper bounds 
above are exact and do not require of confidence in- 
tervals. A method for approximating failure probabil- 
ities that uses both the failure bounds and sampling is 
presented next. 

4.3 Hybrid Method (HM) 

Let be the probability bound associated to the ho- 
mothetic set Ti u (a). In this context 

P[P\ = i>- P [g(U-\u),d) > 0|K(«)] • (15) 

The conditional probability will be estimated via sam- 
pling. Equation 15 requires the generation of samples 
outside 7 T u (a). Since this practice avoids sampling 
the non-failure region of the u- space inside the ho- 
mothet, the accuracy of the approximation is propor- 
tional to the reliability index. If n s is the number of 
samples, an approximation to Equation 15 is 

= v , E": 1 ^b(^- 1 K).d)>0] 

where J, the indicator function, is equal to one if 
its argument is true and zero if it is false, cUj > 0 is 



a weighting factor and u t e 7 i c u (a) is a conditional 
sample. 

A rejection-based algorithm for the generation of 
these samples is presented next. This algorithm gen- 
erates Monte Carlo samples in standard normal space 
and rejects those within the homothet. Its implemen- 
tation is as follows. Instantiate the counter i to one. 

1 . Generate a sample of the standard normal space, 
namely u. 

2. Hyper-spheres : regard this sample acceptable 
if ||w — u|| > Ps- Hyper-rectangles : regard the 
sample acceptable if ||Q(w)|| >\\Q(0n)l 

3. If ft is acceptable, make u, — u and increase the 
counter i by one. 

4. If i < n s + 1 go to Step 1 and repeat. 

An approximation to P\m results from using the 
samples Ui and assuming ay = 0 for all i, in Equation 
16. The number of non- acceptable samples generated 
is proportional to the size of 7 i u (a). Therefore, the 
efficiency of this algorithm diminishes with the value 
of the RI, hence with the value of P[iF]. An alterna- 
tive method, well-suited for cases in which P\m is 
small, was developed but omitted here due to space 
limitations. 

5 EXAMPLE 

A two-dimensional problem in d and p has been 
selected to allow for visualization. The same ex- 
ample is considered in (Crespo, Giesy, and Kenny 
2006), where additional aspects of the methodology 
are exercised. The constraint set is g = [3 d 2 — 4 p\ — 
Ad ^.2 sin (p 2 di - pf), - sin {p\p 2 - sin(2p 1 - 2)) - 
didiPi +p 2 , d 1 +pid 2 -4pip 1 -4:sm(2p 1 -2p 2 ), 
2(Pi + p 2 )sm(pl - d 2 ) - 2p l p 2 {d 2 + 2 p\ - 2) + 
d\ — 6 Pi\ T • The designated point is p = [1, 1] T . The 
boundary of the FDS will be shown in subsequent fig- 
ures as a solid line. In this example, soft constraints 
are considered. In what follows, u = 0 and f p (p) = 
/ P i(pl)/p2(p2), where f p i(pl) = 77(1,0.05) and 
f p2 (p2) = iV(l,0.2). Since the support of this joint 
probability density function is not a subset of largest 
homothetic set found in (Crespo, Giesy, and Kenny 
2006), the usage of soft constraints is justified. This 
implies that the RDS is empty and that P[T \ ' s non- 
zero for all designs. The distribution of — log(P[7 r ]) 
for n s = 10000 samples, generated via Monte Carlo 
Sampling (MCS), is shown in Figure 1. Figures show- 
ing — log(P[7 7 ]) are used to illustrate the behavior of 
the approximation when the failure probabilities are 
small. The “hole” at the interior of the FDS results 
from using a finite number of samples. In the vicin- 
ity of this region, the approximation to P[P] is in- 
accurate and very noisy. For a given number of sam- 
ples, sampling-based RBDO algorithms are not only 
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unable to discriminate among designs attaining small 
failure probabilities, but more importantly, wrongly 
identify zero-failure probability solutions. Besides, 
reliability assessments resulting from sampling lead 
to piece- wise constant functions at every design point 
no matter what n is. Spherical (not shown) and 
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Figure 1: — log(P[jF]) via MCS for n s = 10000. 


substantial advantages would result from using the 
bound ip as an approximation to P[P\. Locally, either 
ips or i J n can be the smallest bound. The First Order 
Reliability Method (FORM) and Second Order Relia- 
bility Method generate approximations that also over- 
come some of the deficiencies of sampling. The HM 
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Figure 3: — log (ips) for hyper-spheres. 
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Figure 2: Rectangular RIs for m = [1, 1] T . 

rectangular (Figure 2) RIs in the FDS were calcu- 
lated. Noticeable differences between the two func- 
tions are apparent. Locally, either the spherical-RI 
or the rectangular-RI can be the largest. Overall, the 
maximum rectangular-RI is greater than the maxi- 
mum spherical-RI. Failure probabilities bounds for 
circular and squared homothets are shown in Figures 
3 and 4. Note that in both cases the bounds are smooth 
functions. Since (i) the calculation of the bound not 
only requires considerably fewer function evaluations 
that MCS (about 50 times less for this example), and 
(ii) the bound is continuous everywhere in the FDS, 
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Figure 4: — log(Lre) for hyper-rectangles. 

is used next. Figures 5 and 6 show — log(P[P]) for 
circular and squared homothets respectively. In both 
cases, n s = 2000 samples are used. The approxima- 
tion not only inherits the nice features of the bound, 
such as its smoothness, but it is also is more accurate 
and efficient than MCS. Overall, less than one fourth 
of the function evaluations used to generate the MCS 
results were used by the HM. As expected, the HM 
approximation compares well to the FORM approx- 
imation when P[J-\ is small. On the other hand, the 
HM outperforms the FORM approximation when the 
failure probabilities are not small. For instance, while 
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Figure5: — log( J P[J ?= '] ) for circular homothets viaHM. 



d, 


Figure 6: — log(P[JF]) for squared homothets via HM. 


the HM leads to P[P\ = 0.623 for d = [-2.2, 0.73] T , 
FORM leads to P[P] = 0.379. The corresponding 
MCS approximation, that for this case can be used as 
a yardstick, is P[P] = 0.614. Therefore, the HM com- 
pares very well to MCS for not-small failure probabil- 
ities values and to FORM for small values. 

6 CONCLUSIONS 

Strategies for analyzing and implementing hard and 
soft constraints in a reliability-based setting are pro- 
posed. Emphasis is given to uncertainty models pre- 
scribed by hyper-spheres and hyper-rectangles. While 
the former geometry can be used to deal with un- 
certain parameters of comparable fidelity, the second 
one allows for the manipulation of parameters with 
different fidelity. The mathematical framework devel- 
oped allows the determination of whether it is feasible 
or not to satisfy hard constraints. Strategies for im- 
plementing hard constraints into design optimization 


schemes are proposed. Probabilistic uncertainty mod- 
els are used in the implementation of soft constraints. 
The methodology permits determination of closed- 
form expressions for upper bounds to the failure prob- 
ability. A hybrid method for the approximation of 
failure probabilities, based on the upper bounds and 
conditional sampling, is proposed. Numerical experi- 
ments show that this method results in substantial im- 
provements in accuracy and efficiency as compared to 
alternative methods. 
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